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ABSTRACT 

We have obtained optical and infrared photometry of the quiescent soft X- 
ray transient XTE Jl 118+480. In addition to optical and J-band variations, we 
present the first observed H- and K^-band ellipsoidal variations for this system. 
We model the variations in all bands simultaneously with the WD98 light curve 
modeling code. The infrared colors of the secondary star in this system are 
consistent with a K7V, while there is evidence for light from the accretion disk in 
the optical. Combining the models with the observed spectral energy distribution 
of the system, the most likely value for the orbital inclination angle is 68° ± 2°. 
This inclination angle corresponds to a primary black hole mass of 8.53±0.60 Mq. 
Based on the derived physical parameters and infrared colors of the system, we 
determine a distance of 1.72±0.10 kpc to XTE J1118+480. 

Subject headings: binaries: close — infrared: stars — stars: individual (XTE 
Jl 118+480) — stars: variables: other — X-rays: binaries 
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1. Introduction 

Transient low mass X-ray binaries (LMXBs) exhibit large and abrupt X-ray and optical 
outbursts that can be separated by decades of quiescence (Chen, Shrader & Livio 1997). 
For these systems, the compact object is a black hole or a neutron star, and the companion 
is normally a low-mass K- or M-type dwarf-like star (Charles & Coe 2003). During their 
periods of quiescence, these systems are faint at X-ray, optical, and infrared (IR) wavelengths. 
While the quiescent X-ray emission can be caused by accretion onto a compact object or 
thermal emission from the surface of a neutron star (Garcia et al. 2001; Wijnands 2004), 
the companion can dominate the luminosity at optical and IR energies. During the binary 
orbit, the changing aspect of the tidally-distorted companion causes a periodic modulation 
of the optical and IR emission (Celino 2001). Measurements of these "ellipsoidal variations" 
provide information about the physical parameters of the binary. 

XTE J1118+480 («2ooo = 11^'18"^10.85% (52ooo = 48°02'12.9") was discovered with the 
all-sky monitor (ASM) on the Rossi X-Ray Timing Explorer by Remillard et al. (2000) on 
2000 March 29, while Garcia et al. (2000) spectroscopically identified its 12.9 magnitude 
optical counterpart. Owing to its position in the Galactic halo, this high latitude {b = 
-1-62°) system has been observed by numerous groups over many wavelength regimes. Recent 
orbital parameters determined from optical spectra suggest an orbital period of 4.078 hr and 
a secondary star radial velocity semi-amphtude of 709±7 km (Torres et al. 2004). These 
values imply a mass function oi f(M)— 6. 3±0. 2 M©, identifying the compact object as a black 
hole. 

Determining a precise black hole mass requires an accurate measurement of the orbital 
inclination angle of the system. As discussed in Gelino, Harrison, & Orosz (2001), the 
best way to find the inclination angle in a non-eclipsing system is to model its infrared 
ellipsoidal light curves. In the IR regime, there is a smaller chance of contamination from 
other sources of light in the system. While modeling several light curves from one wavelength 
regime helps to constrain model parameters, simultaneously modeling light curves that span 
more than one wavelength regime provides tighter constraints than modeling IR light curves 
alone. Previous inclination estimates for XTE J1118-I-480 have come from modeling optical 
ellipsoidal variations as the system approached quiescence (McChntock et al. 2001; Wagner 
et al. 2001; Zurita et al. 2002a). These inchnation angles range from 55° (McChntock et al. 
2001) to 83° (Wagner et al. 2001), and correspond to primary masses of 10 Mq and 6.0 
Mq, respectively. Since this system has been known to exhibit optical superhumps from the 
precession of an eccentric accretion disk on its way to quiescence (Zurita et al. 2002a), it 
is important to determine the orbital inclination angle while XTE Jl 118-1-480 is in a truly 
quiescent state. 
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In order to determine an accurate orbital inclination angle for XTE J1118+480, we have 
obtained B-, V-, R-, J-, H-, and K^-band light curves of the system while in quiescence, 
and simultaneously model them here with the WD98 light curve modeling code (Wilson 
1998). To date, this is the most comprehensive ellipsoidal variation data set published for 
this system. The modeled inclination angle is combined with recently published orbital 
parameters to determine a highly constrained mass of the black hole in this X-ray binary. 

2. Observations & Data Reduction 

We observed the XTE J1118-I-480 field in the optical and IR wavelength regimes. Table 1 
summarizes our observations, while we describe them in detail below. 

2.1. Optical Observations 

We obtained optical observations of XTE J1118-I-480 over three nights in 2003 and four 

nights in 2004 using standard Johnson B, V, and R filters with the 1.5 m telescope at 
the TUBITAK National Observatory (TUG)^ in Antalya, Turkey. The data were obtained 
with the imaging CCD^ on 2003 June 4-6 and the ANDOR CCD^ on 2004 March 18-19 
and 2004 April 23-24. A total of 31, 84, and 83 images were obtained in the B-, V- and 
i?-bands, respectively. The data were reduced with bias frames, dark current frames, and 
dome flat flelds using standard procedures. The instrumental magnitudes of XTE J1118-I-480 
were derived through the use of the PSF fitting algorithms, DAOPHOT (Stetson 1987) and 
ALLSTAR in the MIDAS software package using fifteen PSF stars. In addition, a reference 
star (q;=11:18:07.10 and (5=-|-48:03:53.2) close to our transient source was used to reduce 
the scintillation effects and derive the calibrated magnitudes for a given filter. The reference 
star and target star errors were added in quadrature to produce the final errors on the 
photometric points. The resulting data were combined into 20 bins per orbit. 



^See http://www.tug.tubitak.gov.tr/index.html7eii 

^See http://www.tug.tubitak.gov.tr/aletler/SDSUjmaging_ccd/ 

^See http://hea.iki.rssi.ru/AZT22/ENG/focal.htm 
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2.2. Infrcired Observations 

Infrared data on XTE J1118+480 were obtained on 2003 January 17 and 18 using 
SQIID (Simultaneous Quad Infrared Imaging Device"^) on the 2.1 meter telescope at Kitt 
Peak National Observatory. We simultaneously obtained measurements in the J-, if-, and 
Ks-hands and kept the number of counts in each exposure in the linear regime of each chip. 
The data were reduced using tasks in the upsqiid package in IRAF^. The images taken at 
one position were subtracted from the images taken at a slightly offset position to remove 
the sky, dark current, and any bias level. We then flat fielded the images using twilight sky 
flats. 

Aperture photometry was performed on XTE J1118-I-480 and five nearby field stars. 
Using the IRAF phot package, a differential light curve for each band was generated. We 
followed standard error propagation rules for calculating the differential magnitude errors 
from the target and reference star instrumental magnitudes. The differential photometric 
results show that over the course of our observations, the comparison stars did not vary more 
than expected from photon statistics. We used photometric images of the Hunt et al. (1998) 
AS-11 and AS-18 ARNICA standard star fields, as well as 2MASS field stars to calibrate the 
data set. As with the optical, the data were combined into 20 bins per orbit. 

3. Results and Discussion 

Figure 1 presents the resulting B-, V-, and i?-band light curves of XTE J1118+480, 
while the final J, H, and Kg differential light curves are presented in Figure 2. Despite the 
expectation of detecting superhumps from the 52 day accretion disk precession period, there 
was no evidence of a superhump period when the data were run through a periodogram. 
These results are consistent with the findings of Shahbaz et al. (2005) whose optical data 
taken on 2003 June suggested that if a superhump modulation existed, it was at the <0.50% 
level. Shahbaz et al. (2005) also observed stochastic variability and fast flares in their 
light curves. However, the observed stochastic variability had the same magnitude as their 
photometric error bars, and despite the significant power in the power density spectrum, 
the flares seen in XTE Jl 118+480 had roughly the same magnitude as the spread in the 



"^See http://www.noao.edu/kpno/sqiid/ 

^IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the As- 
sociation of Universities for Research in Astronomy, Inc., under cooperative agreement with the National 
Science Foundation. 
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comparison star measurements. The 2003 and 2004 V- and i?-band light curves presented 
here were consistent in shape and amphtude. As Table 1 shows, two thirds of the nights that 
XTE J1118+480 was observed, data was gathered for at least 75% of an orbit, with more 
than one full orbit covered on four of the nine nights. No evidence for unequal maxima, 
flares, or light curve distortions were found in the unphascd data. The periodogram results 
are consistent with the results from Torres et al. (2004) and therefore all data presented 
here have been phased to their ephemeris. These are the first H- and X^-band detections of 
ellipsoidal variations from this X-ray binary system. 



3.1. Ellipsoidal Models 

The spectral type of the secondary star can be estimated by comparing its red optical 
spectrum with the spectra of stars with various spectral types from the same luminosity class. 
Alternatively, one can use a spectral energy distribution (SED), and published limits on the 
spectral type to not only derive an effective temperature of the secondary star, but to also 
estimate both the visual extinction and contamination level. Given that photometric data 
usually has a higher S/N than spectroscopic data, an effective temperature derived using 
photometry can be just as useful as a spectral type derived from a spectroscopic data set. 
To this end, we compared the observed optical/lR (BVRJHKs) SED for XTE Jl 118+480 
with the observed SEDs for KOV - M4V stars (Bessell & Brett 1988; Bessell 1991; Mikami 
& Heck 1982). We present the best fitting SED in Figure 3, and find that it predicts a 
visual extinction of Ay=0.065±0.020 mags, and a secondary star spectral type of K7V. It 
also includes 60 - 70% hght from the accretion disk at B and 30 - 35% at V. A K5V 
gives a slightly worse fit, and predicts more i?-band light than is observed, as well as a 
smaller amount of disk light in the B- and y-bands, inconsistent with previously published 
values. The visual extinction found here is consistent with the column density adopted by 
McClintock, Narayan, & Rybicki (2004, Nh = 1.2 x 10^'^ cm~^) based on three independent 
measurements, and the spectral type found here is consistent with those found through 
spectral fitting (McClintock et al. 2003; Torres et al. 2004, K5/7 V). 

Light from the primary object or accretion disk in an X-ray binary will act to dihitc the 
amplitude of the ellipsoidal variations of the secondary star. Torres et al. (2004) estimate 
that the secondary star in the XTE Jl 118+480 system contributes roughly 55% of the total 
flux between 5800A and 6400A; however, at R = 18.6, the system was not in true quiescence 
when their observations were taken. Zurita et al. (2002b) determined that XTE J1118+480 
has a true quiescent i?-band magnitude of 18.9. The optical data presented here (i?=19.00) 
are consistent with this value, and thus we contend that the system was in a quiescent state 
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during our observations. In addition, Doppler imaging of the system did not detect any Ho; 
emission from a hot spot or accretion stream (Torres et aL 2004). Since our observations 
took place while XTE Jl 118+480 was in true quiescence, the optical data presented here are 
consistent with both of the Torres et al. (2004) results. 

The most difficult contamination source to extract in the case of XTE Jl 118+480 is the 
one that has the shallowest spectral slope. If the disk contamination in the infrared is based 
on the assumption that the optically thin disk radiates through free-free emission processes, 
and we therefore ascribe the entire V^-band excess to free-free emission, then in the Xg-band, 
the contamination would be ~ 8%. An 8% contamination in the infrared bands would cause 
the observed orbital inclination angle of the system to be underestimated by 2°. What if we 
instead assume that any IR contamination is ascribed to a jet or some form of an advection 
dominated accretion flow (ADAF)? Models by Yuan & Cui (2005) for XTE J1118+480 in 
quiescence show that both the jet and ADAF flux in the IR are predicted to constitute 
significantly less than 8% of the companion flux. 

Irradiation of the secondary star by the accretion disk will affect the symmetry of the 
ellipsoidal light curves. Torres et al. (2004) investigated the possibihty of X-ray irradiation 
powering the Ha emission seen in their Doppler tomograms of the XTE Jl 118+480 system. 
They concluded that it was improbable that this could be the case, and found that the 
strength of the Ha emission from the secondary star was comparable with other K dwarfs. 
Similarly, we sec no significant evidence for irradiation effects in the unphased data or the 
light curves presented here. 

As in Gelino, Harrison, & Orosz (2001), we simultaneously modeled the optical and 
infrared light curves of XTE J1118+480 with WD98 (Wilson 1998). See Gehno, Harrison, 
& McNamara (2001) for references and a basic description of the code, and Gelino (2001) 
for a more comprehensive description. 

The data were run through WD98's DC program which uses a Simplex algorithm for 
initial parameter searches and a damped least squares algorithm for error minimization 
between the data and model light curves. We ran the code for a semi-detached binary with 
the primary component's gravitational potential set so that all of its mass is concentrated 
at a point. The most important wavelength- independent input values to WD98 are listed 
in Table 2. The models were run for a range of inclination angles with parameters for 
a K3V through an MIV secondary star. The secondary star atmosphere was determined 
from solar-metaUicity Kurucz models. We used normal, non-irradiated, square-root limb 
darkening coefficients (van Hamme 1993). We also adopted gravity darkening exponents 
found by Claret (2000), and mass ratios ranging from q = 0.030 - 0.056 (Harlaftis et al. 
1999; Orosz 2001). We assumed that the secondary star was not exhibiting any star spots 
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during our observations. The models were run with varying amounts of additional light in all 
bands. Solving for six consistent light curve solutions simultaneously allowed the rejection 
of many disk light scenarios. 

With 118 degrees of freedom, the best fit model had a reduced of 1.65. While the 
error on the best fit inclination angle was dependent on combining the uncertainties from 
varying all of the model parameters simultaneously, we found that changing the spectral 
type of the secondary from a K5V to an MIV resulted in a change in the orbital inclination 
angle of 1°. Similarly, varying q from 0.030 to 0.056 affected i hj < 1°. We find that the 
best fitting B-, V-, R-, J-, H-, and Ks-hand model has i = 68°, and the parameters found 
in Table 3. Figure 1 presents this model for the optical bands, while Figure 2 presents this 
model for the J-, H-, and Kg-hand light curves. 

3.2. The Model and its Uncertainties 

Based on our optical/IR SED, as well as published values, we adopt a secondary star 
spectral type of K7 (Torres ct al. 2004; McClintock et al. 2001) with a temperature of T^s 
= 4250 K (Gray 1992). The corresponding gravity darkening exponent is f3i = 0.34. 

Modeling six light curves simultaneously is a robust method for constraining the amount 
of extra light in the system. If we model the light curves individually and assume that all of 
the light in the system comes from the secondary star, the best fitting orbital inclination angle 
for XTE J1118+480 varies with wavelength. Since the contaminating light does not have 
the same spectrum as the secondary star the light curves at each wavelength are diluted by 
different amounts, causing the best fit inclination for the affected bands to be different from 
the others. This is the case for the B- and V-band light curves. If a jet or other contaminant 
were to have a fiat spectrum that only affected the IR, the results from the SED would most 
likely not match those obtained from spectral measurements, and we would not be able to 
determine a reasonable parameter set and inclination angle that is consistent throughout the 
data. When fit simultaneously, the best fit inclination angle is i = 68° with 62% disk light in 
B, and 31% disk light in V. These disk light contributions are consistent with those found 
through the SED fitting. 

Using an estimate of 8% for the infrared accretion disk contamination in the system 
gives an inclination of 68^2 '^°, however, if we artificially add in 8% or more infrared light, 
we are unable to obtain a reasonable solution for the orbital inclination angle that is con- 
sistent throughout the entire data set. Therefore, based on the IR colors of the system 
{J-K—l.l±0.3), the SED fit, and the results of the simultaneous light curve modeling, it 
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appears unlikely that the infrared light curves are significantly affected by any such contam- 
ination. In order to determine the final error on the orbital inclination angle, we plotted the 
values as a function of i. Therefore, based on the error in each of the model parameters in- 
cluding q, the spectral type (i.e. temperature) of the secondary star, the amount of observed 
disk light, as well as the photometric error bars, the orbital inclination angle is 68° ±2°. 
We combined the determined inclination angle with the orbital period (P = 0.1699167 ± 
1.72x10-5 d), radial velocity of the secondary star (Kg = 709 ± 7 km; Torres et al. (2004)), 
and and the mass ratio {q — 0.0435 ± 0.0100) to find the mass of the primary object. A 
Monte Carlo routine was used to propagate the errors on the above quantities and gives a 
primary mass of 8.53±0.60 Mq, confirming it as a black hole. 

The constraints on the mass of the black hole in this system presented here represent a 
considerable improvement over those previously published. Wagner et al. (2001) determined 
a mass of 6.0 - 7.7 M© from data obtained before XTE 1118-1-480 had entered a quiescent 
state {R ~18.3). McClintock et al. (2001) gave an upper limit of Mi < 10 Mq, and more 
recently, McClintock, Narayan, & Rybicki (2004) adopted a mass of ~ 8 Mq for their 
calculations of the thermal emission from the black hole in the system. The 8.53 Mq black 
hole mass determined here is consistent with the adopted mass of McClintock, Narayan, 
& Rybicki (2004), and falls nicely into the current observed black hole mass distribution. 
Theoretical models by Fryer & Kalogera (2001) predict that there should exist a greater 
number of 3 - 5 M© black holes than 5-12 Mq black holes; however, most of the determined 
black hole masses thus far have fallen into a 6 - 14 Mq range. In fact, thus far, GRO 
J0422-I-32 is the only system with a compact object that falls into the 3-5 Mq range 
(Gelino & Harrison 2003, 3.97±0.95 Mq). 

Using the mass of the compact object and the orbital period, we computed the orbital 
separation of the two components in the system. We then combined the separation with the 
mass ratio to find the size of the Roche lobe for the secondary star. The temperature of 
the secondary and its Roche lobe radius were then used to find the secondary's bolometric 
luminosity and bolometric absolute magnitude. After accounting for the bolometric correc- 
tion (Bessell 1991), the distance modulus for the J, H, and Kg bands were used to find 
an average distance of 1.72±0.10 kpc. Consistent with results from the modeling and SED 
fitting, this calculation assumes that all of the IR fight in the system originates from the 
secondary star. Table 3 lists all of the derived parameters for the XTE J1118-I-480 system. 
Both the mass and radius of the secondary star are smaller than that of a K7V ZAMS star. 
In addition, the infrared colors of the system and its position in the Galactic halo, support 
the notion that the secondary star in the XTE J1118-I-480 system may be evolved. 
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4. Summciry 

In this paper, we have presented the first observed H- and iC-band elhpsoidal variations, 
as well as the first S-, V-^ and i?-band ellipsoidal variations observed while XTE J1118+480 
was in a truly quiescent state. The derived parameters in Table 3 are based on the modeling 
of these variations including disk contributions of 62% in the B-band and 31% in the V- 
band. Consistent with Shahbaz et al. (2005) and Fitzgerald & Orosz (2003), we do not see 
evidence for any optical superhump light or irradiation in the system. This supports the the 
fact that the system was in a truly quiescent state during our observations. 

While the orbital inclination angle found here is lower than that found by groups who op- 
tically studied the system while approaching quiescence and exhibiting superhumps (Zurita 
et al. 2002a, - 82°), it is consistent with that found from data taken in true quies- 

cence that suggest no significant superhump activity (Fitzgerald & Orosz 2003, i=63 - 73°). 
Furthermore, the distance we find is consistent with those found previously through both 
optical (Wagner et al. 2001; McClintock et al. 2001, 1.9±0.4 kpc, 1.8±0.6 kpc, respectively) 
and infrared (Mikolajewska et al. 2005, 1.4±0.2 kpc) observations. 

Simultaneously modeling multi- wavelength light curves allows us to better constrain the 
amount of disk light in an X-ray binary system. As a result, we have been able to constrain 
the mass of the black hole in the XTE 1118-1-480 system to 8.53±0.60 Mq, and the distance 
to the system of 1.72±0.10 kpc. 

SB acknowledges support from the Turkish Academy of Sciences TUBA-GEBIP (Dis- 
tinguished Young Scientist) Fellowship. EK acknowledges partial support from TUBITAK. 
UK and SB acknowledge the Turkish National Observatory of TUBITAK for providing the 
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Two Micron All Sky Survey, which is a joint project of the University of Massachusetts and 
the Infrared Processing and Analysis Center/California Institute of Technology, funded by 
the National Aeronautics and Space Administration and the National Science Foundation. 
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Table 1. Observations of XTE J1118+480 



Date 


Filter 


Exp Time^ 


% Orbital Coverage 


2003 Jan 17 


J,H,K, 


200 


179 


2003 Jan 18 




200 


172 


2003 Jun 4 


V,R 


240 


14^ ,9^= 


2003 Jun 5 


V,R 


240 


62*^ ,73^ 


2003 Jun 6 


V,R 


240 


44b 40C 


2004 Mar 18 


V 


300 


192 


2004 Mar 19 


R 


300 


224 


2004 Apr 23 


B 


420 


28 


2004 Apr 24 


B 


420 


83 



^Effective exposure time per image in seconds 
^Orbital coverage percentage for the V filter 
'^Orbital coverage percentage for the R filter 
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Table 2. Wavelength- Independent WD98 Input Parameters 

Parameter Value 

Orbital Period (days) 0.1699339 
Ephemeris (HJD phase 0.0)^ 2451880.1086 
Orbital Eccentricity 0.0 

Temperature of Secondary (K) 4250 
Mass Ratio (M2/M1) 0.0435 

Atmosphere Model Kurucz (log g = 4.59) 
Limb Darkening Law Square-root 
Secondary Star Gravity Darkening Exponent /3i=0.34 
Secondary Star Bolometric Albedo 0.676 



^Prom Torres et al. (2004) 
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Table 3. Derived Parameters for XTE J1118+480 



Parameter 


Value^ 


Amount of Disk Light at B (%) 


62±3 


Amount of Disk Light at V {%) 


31±3 


Orbital Inchnation Angle (°) 


68±2 


Primary Object Mass Mi (M©) 


8.53 ± 0.60 


Secondary Star Mass M2 (Mq) 


0.37 ±0.03 


Orbital Separation a (Rq) 


2.67 ±0.06 


Secondary Star Radius Rl^ (Rq) 


0.43 ±0.01 


Distance (kpc) 


1.72 ±0.10 



^Errors are la (Ax^ = 1) 
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Fig. 1. — XTE J1118+480 B-, V-, and i?-band light curves from 2003 (open circles) and 
2004 (filled triangles). The data are plotted over two phase cycles for clarity. Error bars are 
1(7. The solid line represents the best fitting {i — 68°) WD98 model as described in the text. 
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Fig. 2 — XTE J1118+480 J-, and K^-band light curves (circles). The data are plotted 
over two phase cycles for clarity. Error bars arc la. The solid line represents the best fitting 
(i = 68°) WD98 model as described in the text. 
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Fig. 3. — XTE Jl 118+480 phase-averaged optical-infrared quiescent SED dereddened by 
Ay — 0.065 mag (filled circles). Error bars are la. The observed data were compared with 
SEDs for KOV - M4V stars with Ay = 0.045 - 0.085 mag. The best fit SED, normalized at 
H, is that of a K7V with 65% extra hght at B and 33% extra light at V (open stars). 



